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Abstract 

Normal inverse Gaussian (NIG) process was introduced by Barndorff-Nielsen (1997) by sub- 
ordinating Brownian motion with drift to an inverse Gaussian process. Increments of NIG 
process are independent and stationary. In this paper, we introduce dependence between 
the increments of NIG process, by subordinating fractional Brownian motion to an inverse 
Gaussian process and call it fractional normal inverse Gaussian (FNIG) process. The ba- 
sic properties of this process are discussed. Its marginal distributions are scale mixtures of 
normal laws, infinitely divisible for the Hurst parameter 1/2 < H < 1 and are heavy tailed. 
First order increments of the process are stationary and possess long-range dependence (LRD) 
property. It is shown that they have persistence of signs LRD property also. A generalization 
of the FNIG process called n-FNIG process is also discussed which allows Hurst parameter H 
in the interval (n — l,n). Possible applications to mathematical finance and hydraulics are 
also pointed out. 

Keywords: Fractional Brownian motion, Fractional normal inverse Gaussian process, infinite 
divisibility, inverse Gaussian distribution, long-range dependence, subordination 



1 Introduction 

In recent years, there has been a considerable interest on the study of subordinated processes, 
since they have wide applications in modeling of financial time series and hydraulic conductiv- 
ity in geophysics. Also, these processes have interesting probabilistic properties (Sato (2001)). 
A subordinator is a one dimensional Levy process with non-decreasing sample paths (a.s.) 
(Applebaum (2004) p. 49). A subordinated process model formulation was originally consid- 
ered by Mandelbrot and Taylor (1967). They considered a one-sided stable law of index a 
(with < a < 1) for the distribution of the subordinator {T(t),t > 0} with subordination 
to a Brownian motion {B(t),t > 0}, which gives a symmetric stable distribution of index 2a 
for {B(T(t)),t > 0}. They interpreted the subordinator {T(t),t > 0} as trading volume (or 
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number of transactions) upto time t. Clark (1973) used lognormal process as a subordina- 
tor. Madan and Seneta (1990) considered gamma process as a subordinator and they called 
the subordinated process {B(T(t)),t > 0} variance gamma (VG) process. Barndorff-Nielsen 
(1997) introduced normal inverse Gaussian (NIG) process which is obtained by subordinating 
Brownian motion with drift to an inverse Gaussian process for modeling of stochastic volatil- 
ity. Heyde (2005) discussed Student Levy process which is obtained by taking inverse gamma 
process as a subordinator. 

The marginals of these subordinated processes are generally heavy tailed and are more peaked 
around the mode. The subordinated processes are useful for modeling of the returns from 
financial assets (stocks). Indeed, log- returns from financial time series often exhibits a signifi- 
cant dependence structure ( Shephard (1995)). Subordinated processes, with subordination to 
Brownian motion have independent increments. So, the above processes are limited in nature 
and may not be suitable for financial data which exhibits a dependence structure. 
Recently, fractional Brownian motion (FBM) is used instead of Brownian motion for subordi- 
nation, as its leads to LRD property, that is, the covariance function 7*. = cov(X(t), X(t + k)) 
of the stationary time series {X(t),t > 0} decay so slowly such that the series sum J2T=i 7(^) 
diverges (Beran (1994)). FBM is a natural generalization of Brownian motion (Mandelbrot 
(1968)). The FBM {B H (t),t > 0} with Hurst parameter H is a centered Gaussian process 
with B H (0) = and with covariance function 

E{B H {t)B H {s)) = ^[t 2H + s 2H -\t-s\ 2H ], t,s>0, (1.1) 

where a 2 = Var(B H (l)). Also, B H (t) ~ N(0, a 2 t 2H ). The FBM has also the following integral 
representation (Mandelbrot (1968)) 

B H (t) - B H (0) = f7 ^ry{ f_J(t - u) H -" 2 - (-uf-V 2 ]dB(u) + J\t - u) H -" 2 dB{u)}. 

(1.2) 

Its first order increments {X(k) = B H (k + 1) — B H (k)}, called fractional Gaussian noise 
are stationary and exhibits LRD property, when 1/2 < H < 1. Mandelbrot (1997) con- 
sidered {B H (0(t)),t > 0}, where {B H (t),t > 0} is an FBM and the activity time process 
{9(t),t > 0} is a multifractal process with non- decreasing paths and stationary increments. 
Recently, Kozubowski et.al. (2006) introduced and studied fractional Laplace motion (FLM) 
by subordinating FBM to a gamma process. They have pointed out the applications of FLM 
to modeling of hydraulic conductivity in geophysics and financial time series. In the same 
spirit, we subordinate the FBM to in inverse Gaussian process to obtain FNIG process and 
study their probabilistic properties. 
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2 FNIG Process 



The density function of an inverse Gaussian distribution with parameters a and b denoted by 
IG(a, b), is given by 

f(x; a, b) = (2n)- 1/2 ax- 3/2 e ab -^ a2x ~ 1+b2x) . (2.1) 

The inverse Gaussian (IG) subordinator with parameters a and (3 is defined by (Applebaum 
(2004), p.51) 

G{t) = inf{s>0;B(s)+/3s = at}, (2.2) 

where a > and /3 £ R. Here, G(t) denotes the first time when Brownian motion with drift 
(3 hits the barrier at. Another interpretation of IG subordinator is that it is a Levy process 
such that the increment G(t + s) — G(s) has inverse Gaussian IG(at, 0) distribution. We call 
{G(t),t > 0} as IG process with parameters a and (3. Barndorff-Nielsen (1997) introduced and 
studied the subordination of Brownian motion with drift to an IG subordinator. Let Z(t) ~ 
lG(5t, y / 'a 2 — (3 2 ), where < \(3\ < a and 5 > 0. Specifically, he studied the process {N(t),t > 
0}, where N(t) = B(Z(t)) +/3Z(t) + jit and /x G R, called NIG process with parameters a, (3, fi 
and 5, denoted by NIG(a, /3, fi, 5). The process {N(t),t > 0} is a Levy process and its 
marginals have heavier tails than Gaussian distribution. Since NIG process is a Levy process 
its increments do not possess LRD property. Our aim here is to introduce dependence between 
the increments by taking FBM instead of Brownian motion. Let {Bff{t),t > 0} be an FBM 
and let {G(t),t > 0} be an IG process with parameters a > and (3 > 0. Also, we assume 
H = 1/2, when (3 = 0. We now define a new process {X(t),t > 0}, by subordinating FBM to 
an IG process: 

{X(t),t > 0} = {B H (G(t)),t > 0}, (2.3) 

where = stands for equality of finite dimensional distributions. We call {X(t),t > 0} fractional 
normal inverse Gaussian (FNIG) process with parameters 9 = (a, (3, a 2 , H). Note that X(t) 
is a variance mixture of normal distribution with mean zero. Consequently, for alH > 0, X(t) 
can be represented as 

X(t) = aG(t) H Z, (2.4) 

where Z is a standard normal variable and is independent of G(t) having an inverse Gaussian 
IG(cni, (3) distribution. 

Remark 2.1. In the representation \2. 3\) . the distribution of G(t) and Bfj{G(t)) determine 
each other, i.e., 

B H {G{t x )) = B H {G(t 2 )) <{=}► G(h) = G(t 2 ). (2.5) 

The proof is as follows. 

BniGih)) = B„(G(t 2 )) for some t x and t 2 
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<=► G^) 2 " = G(t 2 ) 2 " 
^ G(*i) = G(t 2 ). 

The step above the last one follows by the uniqueness of the Laplace transform. 
Remark 2.2. Since B H (t) ~ N(0, a 2 t 2H ), we have 



E 



exp I iu 



7^ — ;j 



exp l — o u 



\2H 



(2.6) 



2 v t 

v4/so, fry the law of large numbers for the subordinator G(t) (see Bertoin (1996), p. 92), 

G(t) a .s. , a 

> EG{1) = — , as i — > oo. 

6 /J/ 

C/sina toe resuft tfitrf, X n X iffEf{X n ) -> £/(X) /or ever?/ bounded continuous function 
f and taking f(x) = exp(—^u 2 x 2H ) , we have from 112.6}) 



2.1 Basic Properties of FNIG Process 

We first obtain the density of one- dimensional distributions. 

Proposition 2.1. Lei > 0} be an FNIG process with parameter vector (a, (3, a 2 , H). 

Then the density function of X(t) is given by 

at f°° 1 1 x 2 a 2 t 2 

fxti ~ x) = ^ ex P^Pt) ^ - s —exp[--(-^^ + — + P 2 v)\ d V> x e R ' ( 2 - 7 ) 

which is clearly symmetric about origin. 

Proof. Since 

POO 

P{X(t)<x) = / P(B H {y)<x)dF G{t) {y), 
Jo 

we have 



fx(t)(x) = I f BH (y)(x)dF G{t) (y) 

1 x 2 s 1 , „ x l,a 2 t 2 



o 



r°° I x 2 1 .1 err 

= / ^^ exp( -2^ ) 7ir^^ 

= ^exp(apt)^ -^exp[~(^ + ^ + (3 2 y)}dy, x G R, (2.8) 
the result follows. □ 
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Remark 2.3. When (3 = and H = 1/2, 

aat 
n{x' 2 + a 2 a 2 t 



and so {X(£),£ > 0} represents a Cauchy process (see Feller 1971, p. 348). 

Let K v (uj) is the modified Bessel function of third kind with index v. An integral representa- 
tion of K v {u) (J0rgensen (1982), p. 170) is 

K u (uo) = - / u > o. (2.9) 



2 _ .. 

Barndorff-Nielsen (1997) considered the NIG(o;, (3, /i, S) distribution with density 



f(x; a, (3, fi, 5) = a(a, /3, /x, 8)q( X - ) 1 Ki{fog( X ^ )}exp(f3x), 

o o 

where 



a(a,P,ii,S) = ■K~ 1 aexp(8y/[a 2 - /3 2 ) - g(x) = a/(1 + x 2 ), (2.10) 
and Ki is modified Bessel function of third kind and index 1. When H = 1/2, we have from 

f Xt ( x ) = ^l exp (a(3t)(x 2 + aHWy^K^-y/x 2 + aH 2 a 2 ) (2.11) 
7r a 

which is same as the distribution of NIG(/3/er, 0, 0, aat). 

We have the following result concerning the moments of FNIG process. 

Proposition 2.2. If {X(t),t > 0} is the FNIG process with parameters a and (3, then 

E\X(t)\* = c 9 y / |a(^)^- 1 / 2 ^ +1 / 2 exp(a/3t)ir 1/2 „^(a/5t), (2.12) 

where c q = \ — crT( J"^ ) 
V 7T 2 

Proof. Using (12741) . 

£|X(£)|" = <7«E|G(*)| 9ir £7|Z| 9 , 

since G(£) and Z are independent. The r-th moment of IG (at, f3) distribution (J0rgensen 
(1982), p. 13) is given by 

/4 = §=#f^r. rem. (2.13) 



Using Ki/ 2 (uj) = K_i/ 2 (lo) = J -co < e w (J0rgensen (1982), p. 170), we now obtain 

EG f = y|«(^)^- 1 /2 t ^+l/2 exp(a ^ ) ^_ 1/2(a/3t) . 



2i 1 + q 

— IY ) for a standard normal variate 

□ 



The result follows by using the fact that E\Z\ q = 
Z. 

Remark 2.4. Note that for large cu, K v (uS) <■ 

which shows that, E\X(t)\ q ~ c q Q qIi t qH , as t — > oo. This implies the increments Y v (t) 
X(t + rj) — X(t) of the FNIG process behave like the increments of FBM, as the lag r\ 



-p-"uj- 1/2 (J0rgensen (1982) p. 171) 



oo. 



Remark 2.5. Let {X(t),t > 0} be the FNIG process with parameter 9 = (a, (3, a 2 , H). Using 
\2.12) , the kurtosis of X(t) is given by, 



EX(t) 4 
(EX(t) 2 ) 2 



-apt 



apt) 



-1/2 



[apt) 



1/2-2H 



A plot of /?2 (t) against t shows that /?2 is a decreasing function in t. This implies that the 
marginals of X(t) (having /3 2 > 3, for all t > 0) have heavier tails and a larger degree of 
peakedness as compared to the Gaussian distribution. Since /3 2 — > 3, as t — > oo, the process 
X(t) becomes Gaussian with an increasing lag. Observe from Figure 2 that the density of the 
standardized (variance unity) random variable X(l), denoted by f(x) have higher degree of 
peakedness and heavier tails than standard normal density denoted by g(x). Also, the density 
of X(l) tends to a normal density for large values of a. 
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Figure 1: Graph of /3 2 (t) versus t 
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Figure 2: Comparison of pdf of standardized random variable X(l) with parameter 6 = 
(3, 1, 1, 0.75) and the pdf of Z ~ N(0, 1). 



2.2 The Covariance Structure 

The following result concerning the covariance function of FNIG process, follows using the 
fact EX(t)X(s) = \[EX{tf + EX(s) 2 - E(X(t) - X(s)) 2 } and using f[2~T2]) . 

Proposition 2.3. The covariance function of FNIG process {X(t),t > 0} for t ^ s, is given 
by 

2 

EX{t)X{s) = ^aQ 2H - 1/2 \t 2H+1 / 2 ex V (af3t)K 1/2 _ 2H (aPt) 

+s 2H+1/2 exp(aj3s)K 1/2 _ 2H (a(3s) 
_|t _ s \^+i/2 expM()t _ s \))K 1/2 _ 2H (aP(\t - s\))] . (2.14) 

Proposition 2.4. For a = [3, EX(t)X(s) ~ EB H {t)B H {s) as t — > oo, s — > oo and |£ — s| — > 
oo. 

Proof. Without loss of generality, consider the case < s < t, where s — > oo with i — s — > oo. 

Let F(s,t) = 7-r — We discuss the following three cases. 

EB H (t)B H (s) 



Casel. Let - -> 1. Then (- — -) -> 0, and from ffTTjl . EB H (t)B H (s) ~ a 2 t 2H . Hence, using 

~ -^-a[t 1 / 2 e a2 ^ 1/2 _ 2H (a 2 t) + S 1 / 2 e a2s ir 1/2 _ M (a 2 S ) 

v( 27r ) 

^ ~ S yH^ _ ^l/2 e a 2 (t- S )E- 1/2 _ 2H (a 2 (t- S ))j 

-> 1. 

Case Suppose that - — > g for some < q < 1, so that ( — - — ) — > (1 — q), as s — > oo. In 
this case, we have 



a 



Hence, 
F(s,t) 



a- 



EB H {t)B H {s) ~ y [t 2H (g 2 " + 1 - (1 - a) 2 ")]. 



■[t 1 / 2 e a2t K 1/2 _ 2 ^(a 2 t) 



7T g 2 " + 1 - (1 - q) 2Hl 
+ q 2H s^e a2 °K 1/2 ^ H (u 2 s) - (1 - a) 2 ^ - s) 1 ^^ K xl ^ 2H {a\t - s))] 



1. 



Case 3. Suppose that - — > 0, so that ( — - — ) — > 1. In this case, 



EB H (t)B H (s) ~ — t 2H [l-(l 



\2Hl „2tj.2H-1 



a z Hst z 



(2.15) 



2 L ~ v ~ f 

An asymptotic expansion for K v (oj) for large u; (see equation A. 9 of J0rgensen (1982), p. 171) 
is 



2!(8cj) 



3!(8cj) 



where \i = 4i/ 2 . Using the above result, we get \i = (1 — 4if) 2 , and 



a 



, (^-1)^-9) 



= — t 2H (l + ; — — + , , 

w v ; 2 L v 8a 2 t 2!(8a 2 t) 2 



+ s 2ff fi _i_ ^~ 1 , (/^- 1)(/^-9) 

1 8a 2 s 2!(8a 2 s) 2 



+ ■■•) 



fi-1 (p- 1)^-9) 



+ 



(7 



.2H 



(1 



8a 2 (t-s) 2!(8a 2 (t-s)) 2 
» ~ 1 , (/i- l)(/x-9) , 
8a 2 t 2!(8a 2 t) 2 j 



+ ■ 



+ V 1 8a 2 s 2!(8a 2 s) 2 J 

_n _ f >2^ n 1 (^- l)(^-9) 

1 t J 1 8a 2 (t-s) 2!(8a 2 (t-s)) 2 



2 tt,2H-1 

a Hst , 



(2.16) 



(2.17) 
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which follows by using binomial series expansion for (1 ) 2H , and then allowing s, t — > oo 

g 

such that > 0. The result now follows from (12. 15[) and (I2.17p . □ 

2.3 FNIG noise 

Let {X(t),t > 0} be an FNIG process. Then the first-order increment process 

{Y v (t),t > 0} = {X(t + V )- X(t),t > 0}, (2.18) 
for a lag 77 > is stationary. 

We have the following result on persistence of signs LRD introduced by Heyde (2002). A 
process {X(t),t > 0} is said to have persistence of signs LRD property if LRD holds for 
{sgnX{t),t> 0} (see Heyde (2002)), where 

+1, if x > 0, 

sgn(x) = { 0, if x = 0, (2.19) 
-1, if x < 0. 

Proposition 2.5. The increment process {Y v (t),t > 0} has persistence of signs LRD property, 
when 1/2 < H < 1. 

Proof. Note that {sgnY v (t),t > 0} has the same finite dimensional distribution as that of 
{sgn(B h (t + rj) — B h (t)) , t > 0}, which has persistence of signs LRD property (see Proposition 
1 of Heyde(2002)). □ 



We call the process 



{Wj,j EN} = {YMJ - 3 e N} (2.20) 



FNIG noise. Next result on the covariance function of FNIG noise directly follows from ( 12.141) 
and Proposition (12.21) . 

Proposition 2.6. The covariance function of the process {Y v (t),t > 0}, defined by A2.18\) . is, 

for s 7^ t, given by 

EY v (t)Y v (s) = ^lo^) 2 "- 1 / 2 \\t-s + r/| 2H+1 / 2 exp(a(3\t - s + ri\)K 1/2 _ 2H (a(3\t -s + rj\) 
v 2tt p l 



+\s-t + r]\ 2H+1/2 exp(a/?|s - t + r]\)Ky 2 _ 2H {a(3\s - 1 + t?|) 

-2(|* - s\) 2H+1 ' 2 exp(a/?(|t - s\))K 1/2 _ 2H (a(3(\t - a |))] , (2.21) 

and for s = t, 

EY v (t)Y v (s) = E[X( V )} 2 = ( T 2 y|a(^) 2 ^ 1 / 2 r ? 2 ^ 1 / 2 exp(a/?r ? )K 1/2 _ 2H (a/?r / ). 
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The next result shows the asymptotic behavior of the covariance function of {Y v (t),t > 0}. 
Proposition 2.7. Let t > be fixed. Then as s — > oo, 

(2.22) 



EY v (t)Y v (s) ~ a 2 (^) 2H V 2 H(2H - l)(s - t) 2 ^ 2 . 
Proof. Using the asymptotic expansion for K v (uj) given in (12.161) . we have, 



EY v (t)Y„(s) 



a ,a s 



(/i-l)(/i-9) 



.(—) 2H ( s - t -n) 2H (l + 
2 K {3> \> U K 8a(3{s-t- V ) ' 2\{8af3{s - 1 - rj)) 2 



+ (s-t + 7]Y H (l + 



-2(s-ty H (l + 



(/i-1) 



(/i-l)(/i-9) + 



ia^s-t + r/) 2!(8a/3(s - * + 77))* 
(/i-1) ( M -l)( M -9) 



+ 



a/3(s-t) 2\(8af3(s - t)f 



T(f)"(« -*)"[(!- 



+ ■ ■ 



(/i -!)(// -9) | 



+(1 + _^)»(1 



a/3(s-t-T]) 2\(8a(3(s - t - 7])) 2 
(/i-1) ^-1)^-9) 



s -f 



ta/3(s-t + 77) 2!(8a/3(s -t + r/)) : 



+ •••) 



-2(1 



(/i-1) | (/x-l)(/i-9) | 



Sa/3(s-t) 2!(8a/3(s-t)) 2 



v^h^H^H -l)(s- t) 2H - 2 [l + 0((s - f)" 1 )]. 

The last equality follows by using binomial series expansion of (l + ^) 2H and (1 
From fT2~23|) . 



(2.23) 



V ^2H 



EY v (t)Y v (s) ~ a 2 (^) 2 V^(2^ - 1)(« - f) 



s-f 

2H-2 



since lim 0((s — £) ) = 0. 

s— >oo 

Corollary 2.1. FNIG noise process has LRD property for 1/2 < H < 1. 

Proof. The covariance function of {Wj,j > 1} can easily be obtained from (12.211) as 

1 {k) = cov(W j ,W j+k ) -- 



s-f 

(2.24) 
□ 



Z7T P 



(r/(A; + l)) M+1 / 2 exp(«/5r/(A; + l))jr 1/2 _ 2H (a^(A; + 1) 

+ (77(A; - l)) 2H+1 / 2 exp(a(3 V (k - l))K 1/2 _ 2H (a(3 V (k - 1) 
~2( V k) 2H+l / 2 exp(aPr l k)K l/2 _ 2H (af3rik) 



for > 1. 
Further, we have 



7 (jfe) „ a 2 (^) 2H r ] 2H H(2H - l)fc 2H " 2 , as fc - oo. 



From (|235)l . we have ^t(A;) = oo, for 1/2 < H < 1. 



(2.25) 
□ 



fc=i 
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Remark 2.6. Observe that when a = (3 and r] = 1, FNIG noise behaves like fractional 
Gaussian noise with an increasing lag k. Also, the covariances (and correlations) are positive 
when H > 1/2 and negative when H < 1/2, as in the case of FBM. 

A stochastic process {Z(t),t > 0} is said to be self-similar with Hurst parameter H e (0, 1) 
if for any a > 0, we have 

(2.26) 



{Z(at),t> 0} = {a H Z{t),t > 0}. 

A stationary process {Y{t),t > 0} is related to a self-similar process {Z(t),t > 0} by the 
Lamperti transformation defined by 



Z(t) = t H Y(ln(t)), Z(0) = 0, t>0, 



[2.27) 



(see Lamperti (1962)). The process {Z(t),t > 0}, is self-similar with exponent H. The process 
{Z(t),t > 0} may or may not have strictly stationary increments. If we apply the Lamperti 
transformation (12.271) to Y v (t) defined in (12.181) . we obtain self-similar process {W(t),t > 0} 
that do not have strictly stationary increments. In fact, its covariance function for s ^ t, is 
given by 



EW(t)W(s) 



(7 



'^f H ^( st ) H \ln{-) +rf +1 / 2 exp(a/3|Zn(-) + V \)K 1/2 _ 2H (a^\ln(-) + v \) 



a 



2vr 13' 



+\H~) - v\ 2H+1/2 exp(aP\ln(-) - rjl)^ 
s s 



t, 



/2-2H(af3\ln(-) - n\) 
s 

2(|/n(-)|) 2 ^ 1 / 2 exp(a/?(|M-)|))K 1/2 „^(a/?(|/n(-)|))". 
s s s J 



(2.28) 



2.4 Infinite Divisibility 

The generalized gamma convolutions (g.g.c.) were introduced by Thorin (1978a). They 
have turned out to be a powerful tool in proving infinite divisibility of various distributions 
(Bondesson (1979)). A distribution on M + is said to be a g.g.c. if it is the weak limit of finite 
convolutions of gamma distributions (see Steutel and Van Harn (2004), p. 349). Bondesson 
(1979) proved that any density on R + of the form 



M 



A', 



/(*) = 0x^11(1 + J2 c 



x > 0, 



(2.29) 



k=l 



with all the parameters strictly positive, and ctjk < 1, are g.g.c. 's; consequently all densi- 
ties (distributions) which are weak limit of densities of the form (12.291) are g.g.c. 's and thus 
infinitely divisible. Also, if X has a density of the form (12.291) . the weak limit possibilities 
included, then the density of X q , q > 1, is also a g.g.c. (Bondesson (1979)) and hence infinitely 
divisible. 
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Proposition 2.8. When \ < H < 1, the marginal distributions of FNIG process are infinitely 
divisible. 

Proof: It is known that inverse Gaussian distributions or more generally generalized inverse 
Gaussian (GIG) distributions are g.g.c.'s (Halgreen (1979)) and hence infinitely divisible. 
Using Bondesson (1979) result, we conclude that the stochastic variance G(t) 2H has a dis- 
tribution which is a g.g.c, when 2H > 1. Since (B H (G(t))\G(t)) ~ N(0, a 2 G{t) 2H ), and 
X(t) = BH{G(t)) is a variance mixture of normal distribution, it is infinitely divisible (Kelker 
(1971)), when 1/2 <H< 1. □ 

Remark 2.7. (i)Kelker (1971) discussed by an example that the variance mixtures of normal 
distributions can be infinitely divisible, even when the stochastic variance is not. However, we 
have not been able to establish infinite divisibility of FNIG process for < H < 1/2. 
(ii) Halgreen (1979) showed that variance mixtures of normal distribution are self-decomposable, 
if the stochastic variance is a g.g.c. So, the marginals of FNIG process are self- decomposable, 
when 1/2<H< 1. 



2.5 Small Deviation Probabilities 

Small deviation probability (in the logarithmic level) for a stochastic process {X(t),t G T} 
deals with the asymptotic behavior of 

log P {sup \X (t)\ < e}, as e — > 0. 

teT 

For the inverse Gaussian subordinator {G(t),t > 0}, the Laplace exponent <3>, defined by 
E[e- xG{t) ] = e~ m{x \ is given by 

$(x) = a[y/(P + 2x-0\. 

We have 

liminf > 0. 

x^oo Inx 

The following results on the small deviation probabilities for FNIG process {X(t), < t < 1}, 
follows immediately by Theorem 2.1 of Linde and Shi (2004). 

(a) The quenched case (conditional probability := P(-\G(t)))\ 

logP u {\\X\\ q < e} x -$( e - 1/f/ ) a.s. as e -> 0. 

(b) The annealed case (product probability): 

logP{\\X\\ q <e}--^(e^ H ), e -> 0, 

f( x ) fi x ) 

where q G [1, ool , fix) x q(x) means, < liminf ——— < lim sup . . < oo and IIXIL is the 

—0 g(x) ~ x ^o g(x) " M9 

L 9 -norm. 
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3 nth-Order Fractional Normal Inverse Gaussian (n- 
FNIG) Process 



A generalization of FBM called nth-order FBM, is discussed by Perrin et.al. (2001). This 
generalization allows the Hurst parameter H in the range (n — l,n). Note that n-FBM 
is obtained by subtracting terms upto the nth-order of the limit expansion of the kernel 
(t — uf 1 ^ 1 ! 2 from (11. 2J) , leading to the definition of extended FBM: 

Definition 3.1. The nth-order fractional Brownian motion (n-FBM) is defined by 

1 



B n H (t) 



\{t _ u f-l/2 _ (_ u )H-l/2 _ {R _ I)(_ u )H-3/2 t _ . 
2 J ~ J -oo L 2 

1, 

2' 

+ f (t-u) H - l ' 2 dB(u) 



On — 1 4-n—l 
_ *^_^_ u) H-n+l/2_L_ 



n - 1 !J 



tLB(«) 



The covariance function of n-FBM is given by 

n— 1 



^„(t, s ) = (-ir^{|t- s r-^(-i; 

3=0 



2H 

J 



where 



C 



1 



r(2F + l)|sm(7riJ) 



Also, the variance of B\ is 



V[2$(f)] = C£ 



2H-1 
n — 1 



|2# 



^|t| 2 ^(say). 



(3.1) 



(3.2) 



(3.3) 



(3.4) 



To allow Hurst parameter H in the range (n — 1, n), we study nth-order FNIG process by 
subordinating inverse Gaussian process with nth-order fractional Brownian motion (n-FBM). 
We define n-FNIG process as 



{X n (t), t>0}± {B n H (G t ), t > 0}. 



(3.5) 



All the properties of marginals of n-FNIG process follow similar to that of FNIG process by 
replacing o by o n . We can find the covariance function of n-FNIG process with the help of 
(E2D s.t. 

EBl{G{t))EBl{G{s)) = E[(-ir^{\G(t)-G(s)\ 2H -J2(-iy ' 



3=0 

x2H-j , n i \jn(+\2H-j 



J 



G(t) 3 G(s) 3 + G(s) 3 G(ty 



(3.6) 
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Also, we have 

EG(t) j G(s) 2H - j 



k=0 

[i 



E 



k=0 



3 
k 

j 
k 



(G(t) - G(s)) k G(s) 2H - k 
(G(t - s)) k G(s) 2H - k 



xK. 



k=0 



J 
k 



exp (aPt)Kk-i/2(atP(\t - s\)) 



2H-k 



. 1/2 {a/3s){\t-s\) 



fc+l/2 2H-fc+l/2 



(3.7) 



The last statement follows by using (12.1 3p . The covariance function of n-FNIG process can 
be put in explicit form by using (13.71) and (12.131) that is given by 



EX n {t)X n {s) 



2 

ra-l 
j=0 



t-s 



2H+1/2 



exp(a(3\t - s\)K 2 B-i/2(aP\t - s\) 



2H 

J 



7T p 



fc=0 



j 
k 



(exp(a/3t)K fc _i/ 2 (a/9(|t - s 



xK 2H -k-i/2(a/3s)(\t - s \)k+y2 s 2H-k + i/2 + exp ( a p s ) Kk _ 1/2 ( a p(\ t _ s 



xK 2H -k-i/2(af3t)(\t-s\) k+1 / 2 t 



2H-k+l/2 



(3. 



For n = 1, the covariance function of n-FNIG process becomes covariance function of FNIG 
process. Similar to Proposition (12.41) we can show that EX n (t)X n (s) ~ EB H (t)Bjj(s) as 
s — ► oo, t — > oo and It — s\ — > oo. 



3.1 nth-Order Fractional Normal Inverse Gaussian Noise 

For a function /, we denote the kth-order increments by 

A, ( * ) /(x) = A, ( *- 1) /(x + i)-A, ( *- 1) /(x), k>l. (3.9) 

Perrin et.al. (2001) showed that nth-order increments {W H {t) = A\ n) B H (t), t > 0} of n-FBM 
are stationary and hence are called nth-order fractional Gaussian noise. Covariance function 
of nth-order fractional Gaussian noise is given by (see Perrin et.al. (2001)) 

G n ( ?n \ 

cov(W H (t),W H (t + r)) = (-1)"-^ " \r + jl\ 2H , (3.10) 

where r > 0. The nth-order increments of n-FNIG process are stationary and the associated 
noise is obtained by taking nth-order increments of the n-FNIG process. If {X n (t),t > 0} is 
an n-FNIG process, then, for any 77 > 

{Y v n (t),t>0}±{A«X n (t),t>0} (3.11) 
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is a stationary process. We call the process 

{Wj, j EN} = {Y;( V (j - E N} (3.12) 

an nth-order fractional inverse Gaussian noise (n-FNIGN). Using (I3.10p . we find the covariance 
function of n-FNIGN process. 

Proposition 3.1. Let {Y™(t),t > 0} be the n-FNIGN process . Then, for any t, s > 0, where 
s 7^ t, we have 



1 a 2 r 
2 { ~ 



xexp[--[ +p y)\dy. 



where r — t — s. 



4 Simulation 

The sample paths of FNIG process can be simulated by subordinating a discretized IG process 
with fractional Brownian motion. We first simulate a discretized IG process on equally spaced 
intervals. Then we simulate Xk = Bjj(G(tk)), which conditionally on the values of G(tk) is 
a second order Gaussian process with the explicit covariance function of FBM evaluated at 
the values of the IG process. Using Choleski decomposition LL T of the covariance function 
of the random vector BH(G(tk))\G(tk), 1 < k < n. Then X = (Xi, X 2 , • • • , X n ) = LZ, where 
Z ~ N n (0, I n ), a multivariate normal vector. The algorithm is given below. 

(i) Choose a finite interval [0, t], say our aim to have simulated values at to = 0, ti = 
t/n, ■ ■ ■ ,t n _i = (n- l)t/n,t n = t. 

(ii) Generate n IG variates (G\, G2, ■ ■ ■ , G n ), where Gi = G(U), by using the algorithm given 
by Michael, Schucany and Hass (see Cont and Tankov (2004), p. 182). 

(iii) Form the matrix T nxn = (R(Gi, Gj)) , where R(t, s) = \(s 2H + t 2H -\t- s\ 2H ). 

V / nxn 

(iv) Using Cholesky decomposition find L such that T nxn = LL T . 

(v) Generate a n-dimensional standard normal vector Z. 

(iv) Compute X = (X 1; X 2 , ■ • • , X n ) = LZ. Then Xi, X 2 , • ■ ■ , X n denotes n-simulated values 
from FNIG process. 
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2 4 6 8 10 



Figure 3: Sample paths of FNIG process for Hurst parameter H = 0.2. 

Figures 3 and 4 give the simulated sample paths of FNIG process for the case a = (3, and for 
the choices a — .1 (thin lines), a = 1 (medium thick lines) and a = 10 (thick lines). Note 
that sample paths for the case H = 0.8 (Figure 4) have less fluctuations, compared to the 
case when H = 0.2 (Figure 3). This is because of the persistence of signs LRD property of 
the FNIG process for 1/2 < if < 1. 

5 Applications 

Recently, subordinated processes have been found useful models in finance and geophysics. 
Most common processes in literature include variance gamma (Madan and Seneta (1990)), 
normal inverse Gaussian (Barndorff-Nielsen (1997)), FATGBM model (Heyde (1999)), mul- 
tifractal model (Mandelbrot et.al. (1997)), Student processes (Heyde and Leonenko (2005)) 
and FLM (Kozubowski et.al. (2006)). Let P(t) denote the value of a stock at time t. Heyde 
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Figure 4: Sample paths of FNIG process for Hurst parameter H = 0.8. 
(1999) considered the model 

p(t) = P(0) exp(fd + aB(T(t))), (5.13) 

where {T(t),t > 0} is an increasing process with LRD property and having stationary in- 
crements and heavy tails. Also, {B(t),t > 0} is independent of {T(t),t > 0}. Then the 
log-return 

is an uncorrelated process having heavy tails. Mandelbrot et.al. (1997) considered the model 

P(t) = P(O)exp(fl H (0(*))), 

where {9(t),t > 0} is a multifractal process with increasing sample paths and stationary 
increments. It is known that NIG distributions are more suitable than Gaussian or generalized 
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hyperbolic distributions (see Barndroff-Nielsen(1997)). Shephard (1995) observed that log- 
returns R(t) often exhibits dependence structure. Also, the autocorrelations of asset returns 
are often insignificant, but this is not true for small time scales data (Cont and Tankov (2004)). 
Also, Heyde (2002) observed that the log-returns for a financial time series have persistence of 
signs LRD property, i.e., the process has the tendency to change signs more rarely than under 
independence or weak dependence. We have shown that FNIG process has LRD property and 
also have persistence of signs LRD property. Also, FNIG process is a non-stationary process 
with stationary increments, and reduces to Gaussian process, for sufficiently large lags, and 
hence could be useful for large-scale models based on fractional Brownian motion and non- 
Gaussian behavior on smaller scales. Processes with these properties were used by Molz 
and Bowman (1993) and Painter (1996) for analyzing hydraulic conductivity measurements. 
Also, Meerschaert et.al. (2004) used FLM to study the hydraulic conductivity measurements. 
Hence, we believe that FNIG process could serve a useful model in the area of finance, biology, 
geophysics and etc., where the dependence structure of increments and heavy tailedness of 
marginals play an important role. 

Acknowledgment The first author wishes to thank Council of Scientific and Industrial Re- 
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